ABSTRACT
With the emergence and global spread of the COVID-19 pandemic, the scientific community worldwide has focused on search for new therapeutic strategies against this disease. One such critical approach is targeting proteins such as helicases that regulate most of the SARS-CoV-2 RNA metabolism. The purpose of the current study was to predict a library of phytochemicals derived from diverse plant families with high binding affinity to SARS-CoV-2 helicase (Nsp13) enzyme. High throughput virtual screening of the Medicinal Plant Database for Drug Design (MPD3) database was performed on SARS-CoV-2 helicase using AutoDock Vina. Nilotinib, with a docking value of -9.6 kcal/mol, was chosen as a reference molecule. A compound (PubChem CID: 110143421, ZINC database ID: ZINC257223845, eMolecules: 43290531) was screened as the best binder (binding energy of -10.2 kcal/mol on average) to the enzyme by using repeated docking runs in the screening process. On inspection, the compound was disclosed to show different binding sites of the triangular pockets collectively formed by Rec1A, Rec2A, and 1B domains and a stalk domain at the base. The molecule is often bound to the ATP binding site (referred to as binding site 2) of the helicase enzyme. The compound was further discovered to fulfill drug-likeness and lead-likeness criteria, have good physicochemical and pharmacokinetics properties, and to be non-toxic. Molecular dynamic simulation analysis of the control/lead compound complexes demonstrated the formation of stable complexes with good intermolecular binding affinity. Lastly, affirmation of the docking simulation studies was accomplished by estimating the binding free energy by MMPB/GBSA technique. Taken together, these findings present further in silco investigation of plant-derived lead compounds to effectively address COVID-19.
Subject(s)
Methyltransferases/antagonists & inhibitors , Methyltransferases/metabolism , RNA Helicases/antagonists & inhibitors , RNA Helicases/metabolism , SARS-CoV-2/enzymology , Viral Nonstructural Proteins/antagonists & inhibitors , Viral Nonstructural Proteins/metabolism , Antiviral Agents/chemistry , Antiviral Agents/metabolism , Antiviral Agents/pharmacokinetics , Antiviral Agents/toxicity , Binding Sites , Biological Availability , Computational Biology/methods , Databases, Chemical , Drug Design , Humans , Hydrogen Bonding , Methyltransferases/chemistry , Molecular Docking Simulation , Molecular Dynamics Simulation , Phytochemicals/chemistry , Phytochemicals/metabolism , Plants, Medicinal/chemistry , Protein Binding , Protein Domains/drug effects , Pyrimidines/chemistry , Pyrimidines/metabolism , Pyrimidines/pharmacokinetics , Pyrimidines/toxicity , RNA Helicases/chemistry , Structure-Activity Relationship , Thermodynamics , Viral Nonstructural Proteins/chemistry , COVID-19 Drug TreatmentABSTRACT
Novel antiviral active molecule 2- [(4,6-diaminopyrimidin-2-yl)sulfanyl]-N-(4-fluoro- phenyl)acetamide has been synthesised and characterized by FT-IR and FT-Raman spectra. The equilibrium geometry, natural bond orbital calculations and vibrational assignments have been carried out using density functional B3LYP method with the 6-311G++(d,p) basis set. The complete vibrational assignments for all the vibrational modes have been supported by normal coordinate analysis, force constants and potential energy distributions. A detailed analysis of the intermolecular interactions has been performed based on the Hirshfeld surfaces. Drug likeness has been carried out based on Lipinski's rule and the absorption, distribution, metabolism, excretion and toxicity of the title molecule has been calculated. Antiviral potency of 2- [(4,6-diaminopyrimidin-2-yl)sulfanyl]-N-(4-fluoro-phenyl) acetamide has been investigated by docking against SARS-CoV-2 protein. The optimized geometry shows near-planarity between the phenyl ring and the pyrimidine ring. Differences in the geometries due to the substitution of the most electronegative fluorine atom and intermolecular contacts due to amino pyrimidine were analyzed. NBO analysis reveals the formation of two strong stable hydrogen bonded N-H···N intermolecular interactions and weak intramolecular interactions C-H···O and N-H···O. The Hirshfeld surfaces and consequently the 2D-fingerprint confirm the nature of intermolecular interactions and their quantitative contributions towards the crystal packing. The red shift in N-H stretching frequency exposed from IR substantiate the formation of N-H···N intermolecular hydrogen bond. Drug likeness and absorption, distribution, metabolism, excretion and toxicity properties analysis gives an idea about the pharmacokinetic properties of the title molecule. The binding energy -8.7 kcal/mol of the nonbonding interaction present a clear view that 2- [(4,6-diaminopyrimidin-2-yl)sulfanyl]-N-(4-fluoro- phenyl) acetamide can irreversibly interact with SARS-CoV-2 protease.